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Abstract 

We study some of the cosmological imprints of pre-inflationary particles. We 
show that each such particle provides a seed for a spherically symmetric cosmic 
defect. The profile of this cosmic defect is fixed and its magnitude is linear in a 
single parameter that is determined by the mass of the pre-inflationary particle. 
We study the CMB and peculiar velocity imprints of this cosmic defect and suggest 
that it could explain some of the large scale cosmological anomalies. 
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1 Introduction 



In this paper we study the cosmological imprints of a pre- inflationary particle (PIP). By 
a PIP we mean a particle that existed at the latest during the end of the pre-inflationary 
period and the beginning of inflation. The reason we focus on PIPs is that this is the 
simplest and most generic prediction concerning the speculative pre-inflationary era, and 
indeed such particles are created in many pre-inflationary scenarios. 

At first sight there seems to be little motivation for this study. The reason is that one of 
the most basic features of cosmic inflation is that it erases all details of the pre-inflationary 
period. It is this feature that makes inflation such a predictive model. The parameter 
that determines how precise this statement is, is X = e~ AN , where A A" = N tot — Nbb 
is the difference between the total number of e-foldings accumulated during inflation and 
the number of e-foldings required to resolve the big-bang puzzles (horizon and flatness 
problems, etc). Roughly speaking, the amount of pre-inflationary imprints in the visible 
universe scales like a polynomial in X, the details of which depend on the pre-inflationary 
period. For example, imprints of PIPs in the visible universe scale as X 3 (since they get 
diluted like l/a 3 (t)). 

In large field models of inflation AA" is typically large and the probability of finding 
pre-inflationary imprints in the visible universe is negligible. However, in small field 
models of inflation AA" can be fairly small and there might be pre-inflationary imprints 
in the visible universe. 

Pre-inflationary imprints, if they exist, are expected to show up in the largest scales in 
the visible universe. Interestingly enough, some observations at the largest scales appear 

to disagree with acdm Q, uM, H, SB, 0, H, S, Q- 

The results are far from being 
conclusive (partially because of cosmic variance). Still, we find them intriguing enough 
to entertain the possibility that AA" is small, and attribute at least some of them to pre- 
inflationary effects in general and PIPs in particular. In fact, the most generic prediction 
of a small AA" is a significant drop in the power spectrum at the largest scales ll| in 
accord with COBE [l| and WMAP 0, [3] (for a recent analysis see js|). 

In section 2 we show how a PIP affects the evolution of the inflaton during infla- 
tion. We show that to leading order this effect is not sensitive to the details of the 
pre-inflationary period, but only to a single parameter, A, determined by the mass of the 
PIP. In the linearized approximation the dependence is proportional to A, which simplifies 
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the calculations throughout the paper. We also calculate the ideal signal to noise ratio 
{S/N) associated with such a particle. By ideal S/N we mean a situation in which we 
have access to the full cosmological data in the visible universe. 

In section 3 we focus on the PIP's imprints on large scale structure. We show that a 
PIP provides a seed for a spherically symmetric cosmic defect (SSCD), in the form of an 
extended region in space with higher or lower mean density. The shape of the SSCD is 
fixed and its magnitude is linear in A. The sign of A determines whether this is an over 
dense or an under dense region. We study in detail the gravitational potential and the 
profile of the SSCD. In section 4 we study the PIP's imprints on the CMB and investigate 
their dependence on the location of the SSCD. We show that the main effect lies in 
the low-/ modes and that there is a competition between the Sachs- Wolfe (SW) and the 
Integrated SW (ISW) effects induced by the SSCD. We calculate the CMB signal to noise 
ratio as a function of the location of the SSCD and find that the SW-ISW competition 
leads to some interesting effects that could help in identifying the SSCD's imprints in the 
CMB sky. We also show that a SSCD located at a particular distance from us provides a 
possible explanation to the WMAP cold spot. 

The peculiar velocity imprints of the SSCD are calculated in section 5. We find that 
the peculiar velocity signal of a relatively nearby SSCD is larger than its CMB signal. 
In section 6 we consider the possibility that the anomalous peculiar velocity reported 
recently by 0, S] is induced by a SSCD. We calculate the CMB S/N associated with a 
SSCD assuming it is responsible for this anomalously large peculiar velocity as a function 
of its location. We find that typically such a SSCD is CMB detectable, when focusing on 
the right observables. We conclude in section 7. 

2 Point-like particle during inflation 



In this section we generalize the discussion of [13j and study the effects a single Pre- 
Inflationary Particle (PIP) has on the evolution of scalar inhomogeneities. 

For simplicity we assume that the PIP is located at the origin of a spatially flat, 
expanding universe 

ds 2 = -dt 2 + a 2 {t)8 ij dx i dx j . (2.1) 
To calculate the scalar inhomogeneities we find it most convenient to use the spatially 
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flat slicing gauge in which the inhomogeneities are parameterized by 

<f>(x, t) = <p(t) + S<f>(x, t), and <% = 0. (2.2) 

In this gauge the equation of motion for the perturbation in the inflaton in the presence 
of a PIP (see Appendix A.l for the derivation) reads 

«j» + 3H «--^V»*#=-^, (2.3) 

where 

, 1 V 

A = m - -— m, (2.4) 

is a dimensionless effective parameter controlling the influence of the pre-inflationary 
particle. We work in units where c = {StxG)" 1 ^ 2 = 1 and ' means derivative with respect 
to 0. 

Expanding in the slow-roll parameters, the leading contribution to A comes from ml. 
The reason is simple: if the mass depends on the inflaton then there is a direct coupling 
between the PIP and the inflaton, while constant mass means that the interaction is 
indirect, via gravity. 

Eq. ( 12. 31) differs from the standard equation of motion for the perturbation 8(f) by 
the addition of a source term. Being linear, its solution is a sum of a homogeneous 
solution and a particular inhomogeneous solution. The homogenous solution of the above 
equation describes the quantum fluctuations of the inflaton field. These fluctuations have 
a gaussian distribution function with a vanishing one-point function and a non-vanishing 
two point function (power spectrum)@ 

S S(k-k% (2.5) 

ZKj k=a{t)H 

where H is the Hubble constant during inflation. It is this power spectrum that is believed 
to provide the seed for most, if not all, the large scale structure in the universe. 

Imposing the natural initial condition for the inhomogeneous solution (i.e. that it 
scales like 1 1 rait) in the limit rait) <C 1/H) we find that the late time inhomogeneous 
solution induced by the PIP reads (for details see [13l |) 



k 3 \V32 
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(2.6) 

k=a(t)H 



The normalization factor in this equation depends on the normalization convention of the Fourier 

d :, k 

(2tt) 3 / 2 ' 



integral. We work with F(x) = J e^ s F(k) 
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As usual the equation is evaluated at horizon crossing. Note that A could depend on k as 
well. Notice further that, unlike in the homogeneous solution, here we are talking about 
a one-point function and not a two-point function. The non-vanishing one-point function 



Thus, the PIP creates a spherically symmetric defect in the 



depends only on k = 
inflaton field. 

Eqs. (12.51) and (12.61) imply that both 5<ppip and 5<p scale like H. Therefore, for 
the PIP's imprints to be noticeable we need the dimensionless parameter A to be large, 
independently of H. To make this statement more precise we compute the signal to noise 
ratio {S/N) in an ideal situation, in which a PIP is located at the origin of the survey 
and we have full access to comoving modes in the range k\ < k < A^. By signal here we 
mean 5(f>pjp and by noise the quantum fluctuation power spectrum. We have 

where the 4n comes from the angular integration. It is natural to take k2 to be the scale 
at which the linearized approximation breaks down, which is roughly (5 Mpc/h) -1 . At 
the other end, k\ is bounded by the Hubble scale. Hence the best we can hope for is 

1 v / ideal ^ 

which implies that for the S/N to be larger than one we should have A > a/2/3. Since 
the dependence of S/N on the size of the survey, 1/fci, is only logarithmic, this conclusion 
is not changed much for realistic surveys with l/k\ of the order of 100 Mpc/h. 

We conclude that the imprints of a PIP in a sky survey are detectable if |A| > 1. To 
actually be able to detect the PIP's imprints we need to know what to look for. This is 
described in the next sections. 

2.1 Applicability of the approximation 

So far we assumed that there is a clear separation between the effect caused by the PIP, 
8<f>pip, and the standard power spectrum caused by the usual quantum fluctuation of the 
inflaton (12.51) . For large A we expect this assumption to break down and to find mixing 
between the two. 

In other words, we have neglected the backreaction of 5<ppip which for large A will alter 
the power spectrum in a significant way. To estimate when the backreaction of 5<fipip is 
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small we compare it to the standard driving force of the inflaton, V , 

-±-V 2 5<p PIP (r) « V, (2.9) 
a{t) 2 

where 50p/p(r) = A|Mog(r) is the Fourier transform of (12. 6ft . This condition should be 
satisfied at horizon crossing when the power spectrum is determined. Thus we get 

AnV' 

A < — — = 10 5 , (2.10) 
ri d 

where we have used the relation H 2 = V/3 and the COBE normalization (V 3 ' 2 /V = 
5.169- 10" 4 ). 

The values of A that appear in the rest of the paper are at most 10 2 which implies 
that neglecting the backreaction of 50p/p is a good approximation. 

3 Large scale structure 

So far we have concluded that a PIP will create a spherically symmetric defect in the 
inflaton profile. This in turn provides the seed for a spherically symmetric cosmic defect 
(SSCD) whose properties we study here. We start with the gravitational potential and 
discuss its behavior in the small and large scale limits (compared to the comoving Hubble 
scale at matter-radiation equality). We then discuss the energy density profile of the 
SSCD. 

We work in the linearized approximation in which the gravitational potential associ- 
ated with the SSCD is 

®sscD{k,z) = ^%(k)T(k)D 1 (z)(l + z), (3.11) 
where $o is the primordial gravitational potential induced by the PIP (see Appendix A. 2) 

Mk) = JI^.Ak) = ±^ (3.12) 



and T(k) is the transfer function. The evolution of the potential with the redshift (denoted 
by z) in ACDM is encoded in D\(z)(l + z), where D\(z) is the growth function defined 
by 

5 f°° \ _i_ z ' 

D^z) = -Q m E(z) J ■^ydz', with E(z) = y/Sl m (l + zf + Q A . (3.13) 
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Figure 1: The profile of the gravitational potential <&sscD{r) for A = 1. The solid red 
line stands for the potential &sscD(r), calculated numerically (the transfer function was 
calculated via CAMB). We fit the potential numerically in the region 10 Mpc/h < r < 500 
Mpc/h and the result $/^(r) is plotted with a dashed green line. 

For an actual detection it is useful to know the real space profile which is obtained by 
a Fourier transform of $(&). At distances much larger than the comoving Hubble scale 
during matter-radiation equality, ~ 110 Mpc, the transfer function is approximately 1, 
which yields^! 



®sscD(r,z) 



-X- 



—D 1 {z){\+z) log(r). 



(3.14) 



20vr V 

We denote by r the distance from the center of the SSCD and measure it in units of 
Mpc/h. Using the COBE normalization (V 3 ^ 2 /V = 5.169 • 10~ 4 ) and setting z = 0we 
find 

®ssc D (r,z = 0) = |A|C dog (r), (3.15) 



where 



C = ±1.09 -10" 



(3.16) 



The plus (minus) sign is obtained when XV' < (>) for which (13.151) describes an over 
(under) dense spherically symmetric region. If m does not depend on the inflaton then 
A = — iy ra/V and XV' < 0. Therefore, a PIP with a constant mass provides the seed for 
an over dense region. For a PIP to induce an under dense region (or a void), its A should 
be dominated by m with m'V' > 0. We shall use this large distances approximation 
throughout the paper, when relevant. 



3 We assume for simplicity that n s = 1 and that A is a constant. 
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At distances shorter than 110 Mpc, the transfer function is more complicated and to 
obtain an accurate profile of the potential one has to rely on numerical methods. To 
this end we used a standalone version of CAMB 14j with cosmological parameters taken 
from the joint results of WMAP+BAO+SNI [15|: H = 70.1 km/sec/Mpc, Q b = 0.0462, 
Q c = 0.233 and = 0.721. We examined the resulting potential in the range 10 Mpc/h 
< r < 500 Mpc/h and found a good fit to the numerical curve. The fitting function is 



$ Jit (r) = A ■ HT 5 a— log(l + -) + d\, (3.17) 
\ o + r c J 

where the numerical values of the parameters are 

a = 1.1394 ± 0.0023, b = 20.27 ± 0.23 Mpc/h, 

c = 75.44 ± 0.36 Mpc/h, d = -0.8392 ± 0.0038. (3.18) 

Figure 1 illustrates how well (13 . 1 7[) approximates the CAMB result in this range. 

The energy density associated with the SSCD is determined by the gravitational po- 
tential in the standard way 

$(r) = — 1 —V 2 <S>(r,z = 0). (3.19) 
47rGp m 

At distances much larger than l/k eq ~ 110 Mpc, eq. (13 . 1 5[) is applicable and gives 

933 

5{r) * A-^. (3.20) 

Therefore, at large distances the SSCD has a profile identical to that of an isothermal 
sphere, and in particular M(r) grows linearly with the radius M(r) ~ r. 

As expected, the SSCD has a stronger effect at large scales than the standard large 
scale structure, which decays faster with the distance. An example is the well-known 
NFW profile [16J where 5(r) oc r~ 3 at large distances. This, together with the fact that 
it is spherically symmetric (as opposed to a typical cosmic web that includes filaments, 
walls and voids) suggests that if indeed SSCDs are present in the visible universe (and if 
| A | is large enough), one should be able to tell them apart from the standard large scale 
structure. 

At shorter distances we can combine (13. 17)) and (13. 19)) to find the density profile. 
Expanding around r = 0, we get 

p(r) = A (0.96 - O.llr + 0.009r 2 + ...) . (3.21) 

We see that for |A| of the order of 1 the core of the SSCD involves non linear effects. 
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4 CMB imprints 



In this section we study the imprints of a SSCD in the CMB. Our main goal is to calculate 
the CMB signal to noise ratio (S/N) associated with a single SSCD as a function of the 
only two free parameters: its magnitude, which depends on A, and distance from the 
observer which we denote by r^. The dependence of the signal on A is linear. However, 
the dependence on ro is quite interesting and it involves some interesting physics. 



In the first subsection we focus on the Sachs- Wolfe (SW) effect [13]. Calculating the 
SW signal associated with the SSCD can be done analytically and it illustrates some key 
features. For example, we show that, as expected from such a SSCD, the main contribution 
to the S/N comes from the low-/ modes. This means that the late integrated SW (ISW) 
effect cannot be neglected. In the second subsection we add the ISW effect associated 
with the SSCD and calculate the combined SW+ISW signal numerically. We find that 
the combined signal has an interesting dependence on r . Its main feature is a dip in the 
signal at around r ~ 4700 Mpc/h. 

4.1 SW signal 

To calculate the anisotropy we work in the following coordinate system (see Figure 2). 
We take the z-axis (not to be confused with the notation for redshift z) to point towards 
the center of the SSCD. The distance between the SSCD and the observer is denoted by 
r and is measured in units of Mpc/h. The last scattering surface is a sphere of radius 
ri ss from the observer . In principle r can be larger than r/ ss = 9750 Mpc/h. Practically, 
however, this case is not very interesting since it yields a small S/N. 

In this coordinate system the gravitational potential associated with the SSCD is 

^sscoir, z) = XCD^il + z) log(|f - f |), (4.22) 

where C = Dl ^ =0 ^ = 1-425 • 10~ 5 and |r — r | = \fr 2 + Tq — 2rr cos 6. The SW temper- 
ature anisotropy associated with the SSCD reads 

ST W $SSCD{ri ss ,Zlss) 1 / 2 , 2 o Q\ f A oo\ 

-jr- = 3 = q^C log{r lss + r - 2r taa r cos 0) , (4.23) 

where we assume that the decoupling occurs during the matter dominated epoch, which 
means that D\(zi ss ){l + z\ sa ) = 1. In (14.231) we ignore the monopole term ($s , 5cd(0)/3) 
which does not yield any interesting effects. 
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Figure 2: The basic setup. The z-axis points towards the center of the SSCD. The distance 
between the SSCD and the observer is denoted by tq. Due to the rotational symmetry 
of the setup we need to specify only one angle defined with respect to the z-axis, which 
we denote by 9. The distance between the observer and the photon along its path we 
denote by r 7 and the last scattering surface is located at a distance ri ss (Iss stands for 
last scattering surface). 



To calculate the S/N we decompose the temperature anisotropy into spherical har- 
monics 

rrp OO m=l 

yW) = EE aimYim(d, 0). (4.24) 

1=1 m=—l 

The fact that the potential profile (I4.22p has azimuthal symmetry means that we get a 
non- vanishing signal only for m = 0. In other words, all the information about the SSCD 
is in the one-point function of the m = modes 

S t = (a l>m=0 ). (4.25) 

Combining (@~23]), (021 and ffl~2B) we find the SW signal to be 

S SW = X ^ V< 21 + 1 ) f 1 dx p i{x) lQg A _ 2yx + y 2^ (426) 

6 J-i 

where y ee r /ri ss and x = cos 8. The full SW signal to noise ratio is 

/ qSW \ 2 

(V =E(^) 2 /^, (4-27) 

^ ' 1=2 

where C\ is determined by the two-point function associated with the standard effects of 
ACDM 

{a lm a} m ,) = 8 lll 8 mml C l . (4.28) 
10 



As usual, the I = 1 mode is ignored since it is mixed with the Doppler effect. We will 
return to this mode in section 5. 

For each I, (14.261) can be computed analytically. For small y we can expand the result 
to find that the leading term scales like 

Sf w oc -Xy l /Vl, (4.29) 

which yields a S/N that decays with / 

{Sfwy/Q ~ X 2 y 2l l. (4.30) 

Thus, as expected, the main contribution to (I4.27P for small y (r <C ri ss ) comes from the 
low-/s. With a bit of work it can be shown that the leading contribution comes from the 
low-/ modes for any y. 

4.2 SW+ISW signal 

The fact that the main contribution to the S/N comes from the low-/ modes means that 
it is not a good approximation to neglect the contribution of the ISW effect to the signal. 
Namely, we have to consider the combined S/N 

/qSW+ISW\ 2 

(^ r -) =J2(Sr + Sl sw ) 2 /C, (4.31) 

^ ' 1=2 

Since the SW and ISW imprints of a large scale structure on the CMB have an opposite 



signature^, one might suspect that the competition between the two could lead to some 
interesting effects. As we show below, this is indeed the case. 
Decomposing the ISW effect associated with the SSCD 

Q). 2 r%rf^ (4J12) 



T Jiss dr 

into spherical harmonics, we find (using (14.241) and (14.251) ) 

S{ sw = (4.33) 

dzD x {z) (1 - /(*)) J dxP t (x) ln(l + y^zf - 2y^z)x), 



6 An over dense region creates a cold spot when the ISW effect dominates and a hot spot when the 
SW effect dominates. 
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where 

f(z) = - 1 + zdD ^ (4 34) 

J[ } D 1 {z) dz ' 1 } 

r 1 {z) is the position of the photon along its trajectory from the surface of last scattering 

to the observer and y 7 (z) = r 7 (z)/r (see Figure 2). 

For quantitative results we have to rely on numerical calculations. There is, however, 

an important qualitative feature of this integral that can be deduced analytically Along 

the position of the photon there is always a region in which r 7 ~ tq. In this region the 

integral over x in (14.331) is of order 1. Therefore, for small y 

Sl sw (x \y, (4.35) 

where again y = r /ri ss . Comparing ( 14.351) with ( 14.291) we see that for small y the leading 
contribution comes from the ISW effectQ 

As we increase y, the ISW effect decays. This happens because in this case a CMB 
photon passes through the SSCD region when the universe is still dominated by cold 
matter. Since in a matter dominated universe the ISW effect vanishes, the ISW contribu- 
tion arises from a small tail of the potential which the photon experiences at later times. 
Therefore at large y the contribution of the ISW decreases with tq. Simultaneously, the 
anisotropy due to the SW effect grows as the SSCD approaches the last scattering surface. 
At a certain point the ISW becomes negligible compared to the SW. Since S sw and S ISW 
have opposite signs, for each / there exists a point, denoted by r h where the two effects 
cancel each other 

Sf W (r = n) + Sf sw (r = n) = 0. (4.36) 

Similar reasoning implies that r; grows with I. 

Numerical calculations of Sf (tq) and S( (tq) show that indeed both decay with I 
and that for any practical purpose it is enough to account for the first 50 multipoles. In 
addition, the numerical calculation shows that indeed r\ grows with I and that for low-/, 
Ti depends fairly mildly on I. For example, in Figure 3 (a) we plot the lowest multipoles 
of the SW+ISW signal versus r and show that 

r 2 = 4400 Mpc/h, r 3 = 4700 Mpc/h and r 4 = 4980 Mpc/h. (4.37) 



y cannot be too small since all along we are using the long range approximation result (|4.22|) which 
is valid for y > 110 /ri ss ~ 10~ 2 . 
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Figure 3: In (a) we plot the signal (Si + S( sw ) in \iK versus tq for A = 1. We demon- 
strate the cancellation of the SW and ISW effects for different multipoles: quadrupole 
(solid red), octupole (dashed green) and / = 4 (dotted blue). For different / the cancellation 
occurs at different rf. r 2 = 4400 Mpc/h, r 3 = 4700 Mpc/h and r 4 = 4980 Mpc/h. Note 
that r\ grows (slowly) with the multipole number. In (b) we plot the critical parameter 
A C r( r o) versus r for different multipole ranges: 2 < / < 50 (solid red), 2 < / < 10 (dashed 
green) and 6 < / < 50 (dotted blue). This illustrates that the main effect comes from 
lowest multipoles. In order to account for higher multipoles we should have large values 
of A cr (r ). 



Since the main contribution to the signal comes from the low-/ modes, the SW-ISW 
cancellation leads to a dip in the S/N around the low-/ cancellation region, i.e. around 
r = 4700 Mpc/h (we will discuss this feature later). 

Our next goal is to pose constrains on the parameters A and Tq. In particular, we wish 
to know for what range of the parameter A at every r the SSCD is observable in a CMB 
experiment. As we discussed earlier, the signal scales with A. Therefore, for every tq there 
is a critical value of A, denoted by A cr (r ), such that for A > A cr (r ) the signal is larger 
than the noise. Namely, for A larger than its critical value, the imprints of the SSCD in 
the CMB are detectable. In Figure 3 (b) we plot the critical value at each location of the 
SSCD. We also study the contribution of the signal in different multipoles. 

For small r the signal is dominated by the ISW effect, which is proportional to r as 
was shown earlier. Thus the critical parameter behaves as 



A cr (r ) oc ri ss /r . 



(4.38) 



As we increase r the ISW starts to deviate from the linear dependence on r and the 
SW effect begins to suppress the ISW effect. As a result there is a minimum to A cr (r ) 
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Figure 4: (a) A hot spot created by a SSCD with A = 1 in the range ro ~ 4200 
Mpc/h — 5000 Mpc/h. The radius of the spot decreases while tq increases. In (b) we 
demonstrate an attempt to relate the SSCD to the WMAP Cold Spot. We plot two 
published temperature measurements of the Cold Spot (in blue pluses) and fit it with a 
profile created by a SSCD of A ~ —95 and located at tq ~ 4700 Mpc/h (solid blue for 
2 < / < 50 and 6 < / < 50 in dashed red). 

at around Tq = 1500 Mpc/h where its value is A cr (r = 1500 Mpc/h) = 4.7. 

When we increase r further the SW-ISW cancellation becomes more effective and 
there is a fairly clear peak at around r$ = 4700 Mpc/h when we account for all the 
information (i.e. for 2 < / < 50 multipoles). Note that the critical parameter is almost 
identical when we keep only the contribution of 2 < / < 10. This illustrates that the 
main signal comes from the low-/ modes. The very low-/ modes, / < 5, are somewhat 
anomalous (for a recent discussion see and so it is natural to calculate the S/N 
excluding these modes. As Figure 3 shows, this leads to two main effects. First, since the 
signal is considerably weaker, A cr (r ) is significantly larger. Second, the peak is shifted 
to around tq = 5700 Mpc/h, which is in the cancellation region of I = 6, I = 7 and I = 8 
(r 6 = 5400 Mpc/h, r 7 = 5600 Mpc/h and r 8 = 5770 Mpc/h respectively). 

The SW-ISW cancellation leads to another interesting effect that can help in detecting 
the SSCD imprints. In the range 4200 Mpc/h < r < 5000 Mpc/h, (2.55 < z < 3.75) 
the temperature profile associated with the SSCD is highly sensitive to its location, r . 
A SSCD located in this range manifests itself as a localized hot or cold spot in the CMB 
temperature map, depending on the sign of A. A rough explanation for this is the follow- 
ing: as we increase ro, more and more low-/ multipoles approach the cancellation region 
vi and thus almost do not contribute to the signal. Therefore, in this range the signal is 
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dominated by some specific multipole / and this multipole number grows with distance. 
As a consequence, the characteristic scale of the temperature profile becomes smaller. In 
other words, the larger tq is, the smaller the spot is. 

Therefore, the radius of such a spot can be used to determine r . This is illustrated 
in Figure 4 (a) where we have plotted the temperature profile 

f {e) = y £sf w +' 8W Y m (6), (4.39) 

1=2 

associated with a SSCD with A = 1 for different values of r in this range. We see that 
throughout this region the profile is dominated by fairly small angles and that it is quite 
sensitive to r . For r = 4200 Mpc/h, the radius of the hot spot is about 32°. For 
r = 4600 Mpc/h, the radius of the hot spot is considerably smaller, about 16°, and the 
amplitude is weaker. For r$ = 4900 Mpc/h the radius of the spot is tiny, 4.5°, and the 
amplitude is so small that it is overshadowed by a cold ring that peaks at about 20°. By 
the time we reach r$ = 5000 Mpc/h, the ring swallows the tiny hole and we have a fairly 
large cold spot of radius 45°. 

The WMAP cold spot 0, 0, H, EJ might be explained in this fashion. The WMAP 



cold spot is a nearly spherically symmetric region with an approximate temperature ST 



73 [iK at ~ 5° 18(| and an average temperature, 5T ~ —20 fiK at angular radius of 



~ 10°, [20j. This fits well with the imprints of a SSCD with A 95 located at r ~ 4700 

Mpc/h (see Figure 4 (b)). In our case the Cold Spot is surrounded by a hot ring with a 
peak at about ~ 30°. This, as well as improved measurements (perhaps from ACT [2l| 
and SPT |22j), might be used to distinguish the SSCD from other possible explanations, 



such localized void 



12 
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24 ] and a cosmic texture [25]. Unfortunately, the fact 



that the magnitude and shape of the hot ring, unlike the cold-spot, is quite sensitive to 
the very low-/ modes (2 < / < 5) makes it less trustable. 

Outside the SW-ISW cancellation region 4200 Mpc/h < r < 5000 Mpc/h, the main 
signal is due to the low-/ modes. Therefore, the SSCD signature at these locations is a 
mild modulation of the temperature anisotropy and has no prominent pattern (see Figure 
5). Such a signature is not easily detectable in the real-space anisotropy map of the CMB 
but can be seen in its decomposition in spherical harmonics. 
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20 40 60 



Figure 5: We show the pattern of ST (9) created by a nearby (solid red) and remote (dashed 
blue) SSCD, both of A = 1. The former creates a diffused hot spot (the ISW effect is 
dominant) and the latter produces a diffused cold spot (the SW effect is dominant). 

5 Bulk flow imprints 

Another interesting large scale effect caused by the SSCD is the peculiar velocity field it 
induces around it. The peculiar velocity is determined by the gravitational potential 



v(z, r) 



H(z) 



-f(z)V$(z,r) 



(5.40) 



3(1 + z) 2 H$tt r 

and it points towards (away from) an over (under) dense region. In this section we 
calculate the bulk flow associated with a SSCD and determine when it dominates the 
rms bulk flow predicted by the ACDM model. The main motivation for this comes from 
analyses such as 0, 0] that reported an anomalously large bulk flow. 
The observed dipole has two contributions 



jyobserved j-^Doppler _|_ j-jGravity 



(5.41) 



both determined by the gravitational potential and can be straightforwardly calculated 
in our setup. The first contribution is due to the velocity of the local bulk relative to the 
last scattering surface. Hence we get 



with 



j-^Doppler _ 
f 



X 



1 /yrpDoppler | 

dxx — — 

T 



n ■ (v bu i k - vi ss (h)). 



(5.42) 



(5.43) 
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In our case we use the gravitational potential generated by the SSCD (14.221) and (15.401) 
to obtain 

, C /5 n M .\ i 

n ■ Vbuik = X— - - L>i(0) -, 

H r lss \3 ) y 



n ■ v lss = -A - ^ - ) 2 (5.44) 



2C (1-yx) 

1 + - 2yx" 

Note that our local motion yields a pure dipole term whereas the motion of the last 
scattering surface induces higher order multipoles as well. We did not account for them in 
the previous section since their contribution is suppressed by a factor of 1/ y/zi ss compared 
to the SW and ISW effects. Similarly the leading contribution to ]j Do PP ler is due to the 
local bulk motion and not the motion of the last scattering surface 

XC* 

D Doppler ^ Lg4 _ / 5 _ 45 ) 

H r 

The second contribution to the bulk flow comes from the gravitational potential the 
CMB experiences. Namely it is the I = 1 modes of the CMB anisotropy, given at leading 
order by 

gravity = gSW + gISW 

In our case Sf w and S[ sw are determined by (14.261) and (I4.33j) respectively. As follows 
from the discussion in the previous section, both S^ w and S{ sw scale like y in the limit 
of small y and so ]j Gramt y is negligible compared to B D °PP ler in this limit. Therefore, at 
small distances (compared to the Hubble scale) we have 

Vobserved = 0.0385-. (5.47) 
^0 

Recall that we are working in units where c = l=Mpc/h. 

Next we consider the rms bulk flow associated with a region of radius R predicted by 



the ACDM model (for a review see 26|]). For large R we are in the linearized region and 
the average bulk flow reads 

0-0183 

Vrms = (5.48) 

For example, for R = 50 Mpc/h (which is relevant for [7|) we get v rms = 110 km/sec. 

Comparing this with (15.471) we find that the peculiar velocity signal of a SSCD is larger 
than the ACDM noise if 

A > 0.475^. (5.49) 
R 
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Note that our estimate is valid for r > 2R and so the signal is larger than the noise only 
for A > 1. 

6 Peculiar velocity and CMB 

In this section we combine the peculiar velocity imprint and the CMB imprints of a single 
SSCD and argue that quite generically the two are correlated. We also point out a possible 
experimental hint for such a correlation. 

The first step is to find the relationship between A and r assuming that a single 
SSCD is responsible for the anomalously large bulk flow. We use the value reported in 
0], Vbuik ~ 407 km/sec, and for small y we use (15.471) to find 

A py (r ) = ±0.0352r , (6.50) 

where as before r$ is measured in units of Mpc/h. At larger distances we have to take 
into account the contribution of D Gramty as well. As in the previous section, we use an 
analytic calculation of the SW contribution and add the ISW contribution numerically. 
As expected, D Gramty has an opposite sign compared to D D °PP ler and so it lowers the 
magnitude of Jj° bserved In Figure 6 (a) we plot A py (r ) (in dashed-dotted purple). 
The linear behavior at small distances is in accord with (16.501) . while at larger distances 
A py (r ) grows faster, as a result of the partial cancellation between ]j D °pp ler anc l jjGravtty _ 

Next, we wish to study the CMB imprints associated with a SSCD with A = A PV (tq). 
In particular, we would like to know for which values of tq the CMB imprints associated 
with the SSCD are detectable, assuming it is responsible for the bulk flow. The answer 
is presented in Figure 6 in two ways. In (a) we plot A cr (r ) in addition to A py (r ). For 
r in which A cr (r ) > \ PV (r ), the SSCD can induce the large peculiar velocity without 
leaving a significant imprint in the CMB. For r such that A cr (r ) < X PV (r ) its CMB 
imprints are detectable. 

Figure 6 (b) presents this in a more quantitative fashion. The CMB S/N is plotted for 
a SSCD with A = A py (r ) in two cases. In one we include the contribution of all modes 
(/ > 2) to the S/N. In the second we consider the contributions of only the more reliable 
modes, those with 6 < I < 18. The reason for the lower bound was discussed in the 



8 This is expected since we know that when the cosmic defect is well outside the horizon (r ^S> ri ss ) its 
main effect is to generate a superhorizon mode. The contribution of a superhorizon mode to £)° bserved 



is known to vanish [27|, |28|. In appendix B we explicitly check that this is indeed the case. 



18 



2< I < 50 CMB 
2<I<10CMB 
6< I < 50 CMB 
PV 




1000 2000 3000 4000 5000 6000 7000 8000 900C 

r Q [Mpc/h] 
(a) 



2< I < 50 

6<l<18 



1000 2000 3000 4000 5000 6000 7000 

r Q [Mpc/h] 
(b) 



Figure 6: (a): The critical parameter A cr (ro) that gives S/N = 1 for different multipole 
ranges (2 < I < 50 solid red, 2 < / < 10 dashed green and 6 < / < 50 dotted blue) 
and the parameter A py (ro) (dashed- dotted purple line) that gives the observed bulk flow 
v buik _ k m / sec _ g ee magnified window for small distances in the upper-left corner 
of the plot, (b): The CMB S/N associated with the SSCD for X PV (r ). The S/N is 
calculated for different multipole ranges (the solid blue line is for the 2 < I < 50 and the 
dashed red one is for 6 < / < 18). 



previous section. The upper bound is due to the fact that the bulk flow direction is close 
to the galactic plane (e.g. (I, b) = (287° ± 9°, 8° ± 6°) in galactic coordinates according 
to ?|). Ignoring modes with I > 18 (which corresponds to angles under 10°) significantly 
reduces the sensitivity to the galactic noise. For this reason we refer to these modes as 
safe modes. 

As we can see, the signal is smaller than the noise if the SSCD is located relatively 
nearby, i.e. ro < 300 Mpc/h for the full signal and ro < 500 Mpc/h for the signal of the 
safe modes. In this case the SSCD can induce the large peculiar velocity without leaving 
detectable imprints in the CMB. The Shapley supercluster, located at around 140 Mpc/h 
(z ~ 0.046) is believed to be responsible for a good fraction of the bulk velocity {29] . I30I ]. 
and it is a possible candidate for a SSCD0. Needless to say, it should be interesting to see 
whether the Shapley supercluster fits the profile discussed in Section 3. If indeed there 
is a SSCD so close to us then it is natural to suspect that there are others in the visible 
universe and look for their imprints. 

At larger distances, we see that S/N > 1 and so the imprints of a SSCD with A = 
A py (r ) in the CMB should be detectable. There is, however, a dip in the S/N of the 



5 We thank A. Dekel for mentioning this possibility to us. 
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safe modes at ro ~ 5600 Mpc/h (see Figure 6 (b)). This dip in the signal is so extreme 
(it goes down all the way to S/N = 2) that such a SSCD is hard to detect. 

Nevertheless we would like to suggest now that perhaps there might already be hints for 
such imprints in the WMAP data. As follows from the discussion in the previous sections, 
the CMB imprint associated with a SSCD located at r ~ 4400 Mpc/h - 4900 Mpc/h 
would be a spherically symmetric hot or cold spot localized roughly in the same direction 
as the bulk flow. Indeed, such a region is visible in the KpO foreground reduction maslq^ 
of WMAP (see Figure 7). Moreover, [3l| argued that this hot spherically symmetric 
region of A8 ~ 16° of the KpO mask is responsible for a good portion of the observed 
correlations at large distances (which in turn are considerably smaller than predicted by 
the ACDM model). 

It is tempting to speculate that this hot spot and anomalous peculiar velocity are 
induced by a single SSCD. Indeed a SSCD with A ~ 150 and r ~ 4500 Mpc/h can 
explain both anomalies. Keeping only the safe modes, 6 < / < 18, we see (figure 6 (b)) 
that the S/N associated with such a SSCD is ~ 4, which means that it is detectable, but 
not necessarily a prominent signature that would pop up in regular analyses of the CMB 
sky. 

As shown in Figure 7, this hot spot, located at (I, b) ~ (244° — 276°, ±13°) in galactic 
coordinates, is slightly away from the direction of the bulk flow as reported in [7J . This 
can be partially explained if we take into account the statistical fluctuations in the pe- 
culiar velocity mentioned in the last section. Using the values v rms ~ 110 km/sec and 
Vbuik ~ 407 km/sec, we get an upper bound of ~ 14.4° on the statistical error of the 
direction measurement of the bulk flow velocity. Together with the reported measuring 
error from [?J this yields a crude total error estimate for the direction of the bulk flow. We 
should mention that in the last decade alone there have been several reports of bulk flow 
measurements with different results. In Figure 7 we show the recent result of Q, which 
incorporates the results from several large scale surveys, a previous result of a combina- 
tion of surveys 34j and the result of 35] , which is based solely on the SMAC survey and 
lies between the first two and agrees much better with the direction of the hot spot. 

We emphasize that current data is not sufficient to reach a definite conclusion about 
the nature of this hot spot. It is not even clear that a cosmological explanation is needed. 



10 The KpO mask is the most aggressive of high intensity masks used to filter out "bright" sources. It is 
based on the K band map and includes a 0.6 degree radius exclusion area around known point sources, 
amounting in total to 23.2% of the sky. 
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Figure 7: A possible correlation between the anomalously large bulk flow and the origin of 
CMB correlations at large scales: we plot the KpO mask of WMAP and superimpose the 
spherically symmetric region singled out by 31| as a major contributor to correlations at 



large angular scales and the direction of the bulk flow according to three reports [ 7| j34j . 135 
allowing a one a error (combination of measurement error reported by 0, 0, 35] and 
statistical error calculated in section 5). 



In particular, it might be explained via the SZ effect associated with electrons in our 
galaxy [3Ji 



33( | . The nice feature of this explanation is that our location in the galaxy 



fixes the location of the hot spot in the right place to agree with experiment. However, 
the effect appears to be too small by roughly two orders of magnitude. We hope that 
near future measurements will provide a better understanding of this hot spot and of the 
bulk flow. 



7 Summary and discussion 

In this paper we studied the cosmological imprints of a pre- inflationary particle (PIP). 
Let us summarize our findings: 

• A PIP provides the seed for a spherically symmetric cosmological defect (SSCD). 
The shape of the SSCD is fixed and its magnitude is linear in the parameter A which 
is determined by the PIP's mass (see (I2.4p ). 

• A SSCD located within a large scale survey is detectable (S/N > 1) if it was seeded 
by a PIP with |A| > 1. This statement depends only logarithmically on the size of 
the large scale survey. 
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• The CMB imprints associated with a SSCD located at a distance ro from us are 
detectable only for |A| > A cr (r ). In the range of current and near future large scale 
surveys, ro < 500 Mpc/h , we find A cr (r ) to be much larger than 1. Hence it should 
be easier to detect a nearby SSCD via its imprint on structure formation. 

• At larger distances A cr (r ) drops and reaches a minimum at r = 1500 Mpc (which 
corresponds to z ~ 0.575), where A cr = 4.7. At even larger distances A cr (r ) starts 
to grow because of partial cancellations between the SW and the ISW contributions. 

• Due to the SW-ISW cancellation the CMB imprint for 4200 Mpc/h < r < 4900 
Mpc/h is a localized spot in the sky. The size of the spot is highly sensitive to r in 
this range. 

We have also entertained the possibility that a SSCD can explain some of the potential 
large scale anomalies: 

• A SSCD with A ~ —95 located at ro ~ 4700 Mpc/h creates a cold spot of the correct 
angular size and amplitude to fit the WMAP Cold Spot. To study the accordance 
better one needs to know the shape of the Cold Spot with better resolution. 

• A SSCD might explain the anomalously large bulk flow reported in 0, 8]. A nearby 
SSCD, r < 300 Mpc/h, can easily induce a large enough bulk flow velocity while 
having negligible CMB imprints (S/N < 1). A SSCD with |A| ~ 170 located at 
4400 Mpc/h < r < 4900 Mpc/h yields a bulk flow of the observed magnitude and 
its CMB imprints are fairly small. 

• Another interesting fact is the possible connection between the location of the hot 
spot singled out in the galactic plane by |31] and the direction of the bulk flow. A 
single SSCD can explain this correlation neatly: If it is located at ro ~ 4500 Mpc/h 
and it has A ~ 150 it will produce both a hot spot of roughly the same angular size 
as that of [31] and a peculiar velocity field of proper magnitude 

Hi. 

We expect future experiments to provide a wealth of new information. In particular, 
wider frequency bands (as in PLANCK) should provide a better view of the galactic 
plane and beyond it. Therefore, more could be said about the spherical region [31] and 
its possible link to the bulk flow, as well as how good the proposed SSCD explains both 



22 



phenomena. In addition, improved large scale surveys will provide a better bulk flow 
measurement by decreasing the root-mean-square peculiar velocity error. This would 
specify the location of the SSCD more precisely. We could also question the nature of 
non-gaussianities in the CMB and try to trace their origin in the bulk-flow-correlated 
coordinate system. 

To complete the discussion, we raise the question: in which pre-inflationary scenarios 
we expect to find a PIP with a large enough |A| to induce noticeable cosmo logical imprints? 
This is a natural question as A is a dimensionless parameter that is much smaller than 1 
for all standard model particles, and hence the cosmological imprints of a single standard 
model PIP is undetectable. 

A particular pre-inflationary scenario in which such a large A is expected was proposed 
recently in 36| (for earlier and closely related work see 37|]). In this scenario, all kinds of 
particles are produced at the pre-inflationary period (either thermally or via a dynamical 
mechanism such as in 38, [39]). Particles which satisfy m V' < have two appealing fea- 



tures compared to those with m V > 0. First, they are extremely effective in resolving 
the overshoot problem associated with small field models of inflation. Second, they are 
light when they are produced and they become heavier as the inflaton rolls down towards 
the slow-roll region. For example, the mass of the particles discussed in {36] grows expo- 
nentially fast with the inflaton. Therefore, such particles are easy to create and by the 
time the inflaton is in the slow-roll region they can easily have a large enough A to seed 
a noticeable spherically symmetric over dense region. 

Another possible scenario is the following: suppose that the pre-inflationary phase 
includes many PIPs, some with m V < and roughly the same number with mV > 
and suppose that the individual |A| of each of these particles is much smaller than 1. Being 
a scalar field, inflaton exchange would cause particles with the same sign of m to attract 
and particles with opposite sign to repeal. This could lead to clumping of particles with 



the same sign of m which effectively give; 



X e ff = J2i ^i; (the sum is over particles with 



the same sign of m ). Though each |Aj| can be much smaller than 1, with enough particle 
clumping we could have \X e ff\ > 1- Note that in this scenario, unlike in the previous one, 
we can get also a detectable under dense SSCD. 



11 We thank N. Afshordi for pointing to us the potential importance of this interaction. 
12 Our current rough estimate indicates that typically the particles do not have enough time to clump. 
This interesting dynamical process deserves however a more detailed study. 
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A Derivation of (E5) and ( f37T2l ) 



In this section we calculate the primordial potential $o(&) which serves as an initial 
condition for the metric perturbation in the post- inflationary epoch. First, we derive the 
modified equation of motion for the inflaton in presence of the massive point-like particle, 
using the spatially flat slicing gauge. Next we evolve the perturbation through the horizon 
and find the primordial potential $o(&) i n the Newtonian comoving gauge. 

A.l The Spatially Flat Slicing Gauge 

Let us start with deriving the modified equation of motion for the inflaton in the spa- 
tially flat slicing gauge. This gauge appears to be one of the most suitable to deal with 
perturbation theory during inflation, since the only approximation that has to be made 
is the slow-roll approximation. The perturbated metric in this gauge reads 

ds 2 = -(1 + Afdt 2 + 2a(t)B \idx' [ dt + a(t) 2 h ij dx i dx j where hij = 8 i5 , (A.51) 

and the perturbed scalar field is 

<f>{x,t) = <f>{t) + S<f>{x,t). (A.52) 

First, we discuss the equation of motion for the inflaton only. In this gauge the 



equation of motion (see 40[ for derivation) reads 



54> + 3H6<P - ~^V 2 5<P — \V" + 2^(3H - ?= + 2^) I 5<p. (A.53) 
a[t) z \ H H 

The R.H.S. of the equation is close to zero in the slow- roll approximation, hence 

54> + 3H5<p - -^- 7 rV 2 54> = 0. (A.54) 
a(t) 2 



The key feature of this gauge is that the equation (1A.54|) does not ignore the coupling 



of the inflaton to the curvature perturbation as it does, for instance, in the Newtonian 
conformal gauge. 
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We will also make use of one of the perturbed Einstein equations, namely 

which couples the metric perturbation to the perturbation in the inflaton. 

If in addition to the inflaton field there exists a PIP, which we assume to be a non- 
dynamical, massive point-like particle, the following term should be added to the action 

S P artide = - J dmi = J dx°^/-g 00 m, (A.56) 

where we have used dr = \/—goodt = (1 + A)dt. To first order in perturbation, the 
variation of S par ticie with respect to the inflaton is 

SS particle = J d 4 x (j^St + mA^j 5 3 (x) = J d 4 x^g- > ( A -57) 

where 

dm IV' , . i 

A = w - - y ,n. (A.58) 

This additional term modifies the equation of motion of the inflaton so that the complete 
equation reads 

.im^^V^-gj. (A.59) 
The solution to this inhomogeneous equation was discussed in detail in [3] and is 

Perturbations generated during inflation are related to the post-inflationary quantities 
via a gauge invariant variable that is conserved on super-horizon scales 

£(k) = S-8(/)(k) - (A.61) 

In the gauge under consideration, the perturbation in the space-space metric vanishes 
\& = 0. Therefore the conserved quantity at the end of inflation reads 

t(k) = -S-Sm- (A.62) 
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A. 2 Primordial Potential 

After inflation ends it is convenient to switch to the Newtonian conformal gauge, which 
in the absence of the anisotropic stress has the following perturbed metric 

ds 2 = -(1 + 2§)dt 2 + a(t) 2 (l - 2$)dx i dx j 5 ij , (A.63) 

where $ is the gravitational potential. In this gauge the primordial perturbation in the 
gravitational potential is 

*o(fc) = -^(k). (A.64) 
For the particular solution ( 1A.60I) . the primordial potential reads 

*»« = (A - 65) 



where the constant 

H 1 V 3 / 2 

= = . — = 1.984 • 1(T 5 (A.66) 

12^ V 1 ; 

comes from the COBE normalization: = 5.169 • 10~ 4 . 

B Superhorizon Mode - A Consistency Check 

As was shown by the authors of j^] for a CDM universe and [28] in the ACDM case, the 
observed dipole induced by any superhorizon perturbation should be very small. Here we 
calculate the dipole created by a SSCD located far beyond the surface of last scattering 
and show that it is indeed negligible. We do this in both cases (with and without the dark 
energy component). In each case we expand the dipole in powers of r; ss /r and calculate 
the linear effect only. 

First, let us neglect the dark energy. In this case the dipole consists of the peculiar 
velocity term and of the gravity term, which in the CDM case is the SW contribution 
only 

j-jObserved j-yDoppler _|_ j^Gravity 67) 



Each term reads 



gravity = gSW = ^^^J^^ (R6g) 
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and 



J^Dopph 



2C hr 



|-A(o)-^ 

■J oniss 



1 



(B.69) 



In (IB.69[) the last term in brackets is the contribution due to the motion of the last 
scattering surface, which we have approximated by 



v (e = o) lss - v 



v lss ■ n 



7T 



. Iss 



Keeping in mind that in a standard CDM scenario r/ s 



ler 



3 V 3 



2 /H and D l (0) 
1 



(B.70) 
1, we have 

(B.71) 



VlTOOj H r 

We see that the SW contribution and the Doppler term due to the peculiar motion of an 
observer (the first term in brackets) cancel each other. Therefore, the leading order to 
the observed dipole in powers of rj as /ro comes solely from the motion of the plasma at 
decoupling. This effect is only ~ 10~ 2 of each of the cancelled components. In other words, 
the dipole produced by a distant SSCD is suppressed and would be barely observable. 

The same conclusion can be made for the ACDM background. We will show that, in 
agreement with 28j, the first order contribution to the dipole vanishes in the presence of 
dark energy 

j-jObserved j-^Doppler _|_ j^Gravity g ^-g 72) 

For the background under consideration, the gravity term consists of the SW and the ISW 

effects D Gravif y = Sf w + S{ sw . The expression for the SW term is identical to flR68j) . In 

order to calculate the ISW contribution we should integrate over the trajectory of a CMB 

photon from the last scattering surface to the observer. For a superhorizon structure, 

each r' on the trajectory of a CMB photon satisfies r' < ro, which allows us to find an 

I L 

analytical expression for the ISW term. Following [28], we find 

1 



ISW 



#oV3 



^-5 + ^(0) + ^ 

6 6 orliss 



[B.73) 



The expression for the Doppler term is identical to (1B.69j) . Adding up all the contributions 
we obtain ]j° bserved = o, i.e. in the ACDM universe the linear contribution to the observed 
dipole due to a superhorizon SSCD vanishes. 

A consequence of the above analysis is that in order to create the anomalous bulk flow 
by a single superhorizon SSCD, the value of A should be enormous. It is hard to imagine 
any physical process in charge of PIPs with such large A. 
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